Evidence for higher order QED in e'^e pair production at RHIC 
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A new lowest order QED calculation for RHIC e'^e" pair production has been carried out with 
a phenomenological treatment of the Coulomb dissociation of the heavy ion nuclei observed in the 
STAR ZDC triggers. The lowest order QED result for the experimental acceptance is nearly two 
standard deviations larger than the STAR data. A corresponding higher order QED calculation is 
consistent with the data. 
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In a recent publication data were presented by the 
RHIC STAR collaboration on "Production of e+e" pairs 
accompanied by nuclear dissociation in ultrapcriphcral 
heavy- ion collisions" 1] . In this publication the authors 
compared the experimental cross section a with a lowest 
order QED calculation agsD and concluded that "At 
a 90% confidence level, higher-order corrections to the 
cross section. Act = a — gqed must be within the range 
-Q.haQED < Act < 0.2aQED'" ■ 

In this Letter, I will present a new lowest order QED 
cross section calculation result significantly larger than 
that presented in the STAR publication. This new low- 
est order QED cross section is now two standard devia- 
tions above the STAR experimental result. Furthermore, 
I find that a higher order QED calculation gives good 
agreement with the experimental result. 

The STAR experimental cross sections are for a final 
state with an e~^e~ pair in combination with mutual nu- 
clear excitation of both ions by the ions' Coulomb fields 
observed in the Zero Degree Calorimeters (ZDC). Thus 
impact parameter dependent probabilities have to be cal- 
culated both for e~^e~ pair production and for mutual 
Coulomb dissociation of the ions. 

The calculations presented here make use of the 
method of calculating the impact parameter depen- 
dent probability of producing e+e~ pairs previously 
described 3. Use of an appropriate gauge allowed the 
electromagnetic fields of colliding ions to be represented 
as colliding delta functions in the ultrarelativisitic limit. 
The solution of the time-dependent Dirac equation for 
lepton pair production was thereby obtained in closed 
form. This solution contains higher order QED effects 
and goes to the known perturbative (lowest order) result 
in that limit. It is well known that photoproduction of 
e~^e~ pairs on a heavy target shows a negative (Coulomb) 
correction proportional to that is well described by the 
Bethe-Maximon theory The calculated higher order 
heavy ion pair production cross sections show analogous 
reductions from perturbation theory. 

The impact parameter (b) dependent amplitude 
presents a particular numerical challenge since it involves 
the a rapidly oscillating phase exp(ik • b) in the integral 
over the transverse momentum k transfered from the ion 
to the lepton pair. The usual method of evaluating the 



perturbative impact parameter dependent probability is 
to first square the amplitude and then integrate over the 
sum and difference of k and k'. Here I have integrated 
before squaring, and I deal with the rapid oscillations 
with the piecewise analytical integration method previ- 
ously described [2], in particular Eq (A5-A7). In that 
previous b dependent calculation of the total cross sec- 
tion, half of the contribution comes from b > 5000 fm 
and contributions up to 6 = 10^ fm are considered. Due 
to the large values of b contributing, that calculation was 
somewhat crude. However integration over b reproduced 
the known cross sections calculated with the b indepen- 
dent method or calculated from the very accurate analyt- 
ical Racah formula[l] to about 3%. It can also be noted 
that the computed perturbative and eikonal 6 dependent 
probabilities in that paper were in relatively good agree- 
ment with calculations in the literature Q available for 
b < 7000 fm. 

In the present calculation although k takes on large 
values scaled to the momenta of the produced pairs, 
the range of of b considered is limited to under 100 fm 
by the mutual Coulomb dissociation requirement of the 
ZDC tagging. The calculated high momentum piece of 
the cross section was taken as described by the STAR 
experimentalistsfH 0|. Electron and positron tracks were 
required to have transverse momentum px > 65 MeV/c 
and pseudorapidity |?7| < 1.15. In addition, at least one 
of the tracks was required to have momentum p < 130 
MeV/c. Pairs were required to have masses 140 MeV 
< Mee < 265 MeV, pair transverse momentum Pt < 100 
MeV/c, and pair rapidity |y| < 1.15. In the numerical 
calculations the mesh size (piece) over which the analyt- 
ical k integrations is taken was chosen such that dou- 
bling the number of mesh points (halving piecewise size) 
changed the computed result by less than 0.1%. 

Several exploratory calculations were carried out to 
check the accuracy of the lowest order e+e" calcula- 
tions in the high momentum ranges here considered. The 
same expression for the pair production amplitude, de- 
rived in the ultrarelativistic limit, is the basis for both 
impact parameter independent 01 and dependent f5] cross 
section expressions. The impact parameter independent 
calculation of the untagged cross section in an approx- 
imate STAR acceptance range carried out by Hencken, 
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Baur, and Trautmann[8|, 0.322 b, was reproduced with 
a niodified version of my impact parameter independent 
programQ, and with a version of the codeLPAIR[i,[l3], 
both to one percent accuracy. 

The nuclear breakup was calculated using a phe- 
nomenological model based on photodissociation data as 
imput for Weizsacker- Williams calculations ll| for mu- 
tual Coulomb dissociation JJj that was found to suc- 
cessfully reproduce cross sections for mutual Coulomb 
dissociation observed in the zero degree calorimeters at 
RHIC 13|. This phenomenological approach to the nu- 
clear breakup was further incorporated in successful cal- 
culations of t agg ed vector meson production probabilities 
at RHIC [U. 

In analogy with the previous tagged vector me- 
son calculations IJ] and STAR equivalent photon 
calculations [H, the cross section here is computed from 
the product of the pair production probability Pee{b), 
the Coulomb dissociation probability Pxx{b), and a fac- 
tor exp[— P„„(&)] to exclude events where hadronic inter- 
actions occur 

CT = 27r / PxMPee{h) cxp [-P„„ (6)] 6d6. (1) 



authors made use of the realistic Coulomb dissociation 
phenomenology ll|, [ij, observed in the ZDC trig- 



P'xxib) is the unitarized probability that both colliding 
nuclei suffer Coulomb dissociation 



Pxx(b) - [l-exp(-Pc(&)]' 



(2) 



with Pc{b) the non-unitarized probability of a single 
Coulomb excitation. The hadronic interaction proba- 
bility Pnn{b) is calculated in the usual Glauber manner 
from impact parameter dependent nuclear density (with 
Woods-Saxon shape R=6.38 fm a=.535 for An) overlap 
and a nuclcon-nucleon cross section of 52 mb for the 
RHIC center of mass energy of 200 GeV. In calculat- 
ing the pair production probability Pee{b) an analytical 
elastic form factor was employed 
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(3) 



with a hard sphere radius r — 6.5 fm and Yukawa poten- 
tial range a = 0.7 fm. This form very closely reproduces 
the Fourier transformation of the Au density with the 
Woods-Saxon parameters mentioned above[16j. Test cal- 
culations with an alternative form factor that has been 
utilized in the literature [l7j| for analytical simplicity, the 
dipole form with A = 83 MeV/c, gave results in agree- 
ment with those making use of the analytical Woods- 
Saxon form to one percent. 

The lowest order QED calculation (1.9 mb) used by 
STAR for comparison with data (1.6 mb) mentioned 
above made use of an approximate analytical treatment 
of the Coulomb dissociation Q combined with lowest or- 
der QED for the pair production. An alternate low- 
est order calculation (2.1 mb) presented by the STAR 



gers such as adopted in the calculations here but with 
an equivalent photon QED calculation of the e~^e~ pairs. 
Both these calculations left something to be desired. The 
lowest order QED calculation, made use of an overly sim- 
plified approximation to the Coulomb dissociation phe- 
nomenology and involved choosing an arbitrary lower 
cutoff for the impact parameter integration. The equiva- 
lent photon calculation, did not include photon virtuality. 
The QED calculations presented here, both lowest order 
and higher order, have neither of these shortcomings. 

It is useful to change variables for the impact parame- 
ter integration in order that the contribution to the cross 
section is a smooth function of a variable that goes to 
zero at the high and low end. Choose a new variable 
u — {bi/b)^^'^ and for convenience set bi equal to 13 fm. 
Figure 1 shows the contribution to the cross section as 
a function of u. It is interesting to note that although 
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FIG. 1: Tagged pair production cross section as a function 
of u. Corresponding values of the impact parameter b are 
indicated next to tick marks at the top. The dashed line 
is for the lowest order QED calculation, the solid line for 
the higher order QED calculation, and the dotted line for an 
eikonal higher order QED calculation. 

the original simple eikonal inclusion of higher order ef- 
fects gives the same result as perturbation theory for 
untagged pair production, the weighting of low impact 
parameter pairs due to the Coulomb dissociation of the 
STAR acceptance lead to an eikonal higher order result 
that differs from perturbation theory. Between & = 46 
and 6 = 60 the eikonal result crosses over from being be- 
low the perturbation theory result at low h to above the 
perturbation theory result at high h. This crossover has 
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been noted at about 3000 fm in numerical calculations of 
the total e+e" pair production probability at RHIC[3,0|- 

The integrated lowest order QED cross section is 2.34 
mb for the STAR acceptance. This is nearly two standard 
deviations higher than the STAR experimental result of 
1.65 ± 0.23(stat)±0.30(syst) mb[i|. The corresponding 
higher order QED calculation result of 1.67 mb is con- 
sistent with the STAR data. Thus the STAR data seem 
to provide some evidence of higher order QED in ultra- 
relativistic heavy ion reactions. The integrated eikonal 
result, 1.80 mb, is also significantly reduced from the low- 
est order perturbative result. 

I would emphasize that the higher order QED result 
here was calculated with the same physical cutoff of the 
transverse potential that restored the Coulomb correc- 
tions to calculations of the total cross section f^, qualita- 
tively consistent with the approximate higher order cal- 
culations carried out by I vanov, Schiller and Serbo[l^ 
and by Lee and Milstein[l^, In particular, the 

higher order QED effect from the physical cutoff was 
shown numerically to go over to the Coulomb correc- 
tion result of Lee and Milstein in their small k expansion 
limit |21|]. The eikonal result makes use of the original 
formulation of Baltz and McLerran[23| and of Segev and 
Wells [23L (first calculated numerically by Hencken, 
Trautmann, and BaurQ) that was not consistent with 
the more proper regularization of Lee and Milstein. 
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FIG. 2: pt spectra for the produced positrons and electrons. 
Squares are the number of positrons and triangles electrons. 
The dashed line is the lowest order QED calculation and the 
solid line higher order QED. 



QED. The solid line represents the corresponding ex- 
pected numbers from the higher order QED calculation. 
Obviously the higher order calculation is in good agree- 
ment with the data, while the lowest order calculation 
overpredicts. 
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FIG. 3: The pair equivalent mass distribution. 

Figure 3 shows pairs cross section data (filled circles) 
as a function of the pairs invariant mass to be more con- 
sistent with the solid line higher order QED calculation 
than with the dashed line lowest order QED. Likewise the 
data shown as a function of the pair momentum in Fig- 
ure 4 show good agreement with the higher order QED 
calculation but are below the lowest order calculation. 
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Figure 2 is a comparison of the px spectra for the 
number of produced positrons and electrons. A total 
of 52 events was observed in the STAR acceptance. The 
dashed line represents the number of expected positrons 
or electrons within each bin calculated in lowest order 



FIG. 4: The pair Pt spectrum. 

In these calculations it has been as crucial to get the 
overall normalization of the lowest order QED calculation 
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correct as to calculate the reduction due to higher order 
effects. This normalization depended on reliable QED 
pair production calculations, reliable Coulomb dissocia- 
tion calculations, and faithful reproduction of the STAR 
experimental momentum cuts. But if these calculations 
are accepted as reliable, then the above comparison of 
the QED calculations with the STAR data provides the 
first evidence of higher order QED in rclativistic heavy 
ion reactions. 

In the planning stages of RHIC a workshop was held 
at Brookhaven on the topic "Can RHIC be used to test 
QED?" ^ A recent paper concluded, "We think that af- 
ter about 17 years the answer to this question is 'no'." [i^ 
The present results indicate that the answer may turn out 
to be 'yes'. 
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